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Abstract 

A short review of theoretical and experimental results on fragmentation in e + e~ 
annihilation is presented. Starting with an introduction of the concept of fragmen- 
tation functions in e + e~ annihilation, aspects of scaling violation, multiplicities, 
small and large x, longitudinal, gluon, light and heavy quark fragmentation are 
summarized. 
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1 Introduction 



Fragmentation functions are dimensionless functions that describe the final-state single- 
particle energy distributions in hard scattering processes. The total e + e~ fragmentation 
function for hadrons of type h in annihilation at cm. energy y/s, via an intermediate 
vector boson V = j/Z°, is defined as 

F h (x, s) = —^-(e + e- -> V -> hX) (1) 

where x = 2Ef l / y/s < 1 is the scaled hadron energy (in practice, the approximation x = 
x p = 2p h /y/s is often used). Since each hadron of type h contributes to the fragmentation 
function ([]]), its integral with respect to x gives the average multiplicity of those hadrons: 

(n h (s))= f 1 dxF\x,s). (2) 
Jo 

Similarly, the integral of xF h (x,s) gives the total scaled energy carried by hadrons of 
that type, which implies the sum rule 

V f 1 dxxF h (x,s) = 2. (3) 
h Jo 

Neglecting contributions suppressed by inverse powers of s, the fragmentation function 
(0) can be represented as a sum of contributions from the different parton types % = 
u,u } d } d, . . . , g: 

F h (x, s) = J2f 1 -Ci(s; z, a s )Df(x/z, s) . (4) 
i Jx z 

where D\ are the parton fragmentation functions. In lowest order the coefficient function 
C g for gluons is zero, while for quarks Cj = gi(s)5(l — z) where gi{s) is the appropriate 
electroweak coupling. In particular, g^s) is proportional to the charge-squared of parton i 
at s Mf , when weak effects can be neglected. In higher orders the coefficient functions 
and parton fragmentation functions are factorization-scheme dependent. 

Parton fragmentation functions are analogous to the parton distributions in deep inelastic 
scattering (Section 9 of In the fragmentation function x represents the fraction of 
a parton's momentum carried by a produced hadron, whereas in a parton distribution 
it represents the fraction of a hadron's momentum carried by a constituent parton. In 
both cases, the simplest parton-model approach would predict a scale-independent x 
distribution. Furthermore we obtain similar violations of this scaling behaviour when 
QCD corrections are taken into account. 



2 Scaling violation 

The evolution of the parton fragmentation function Di(x,t) with increasing scale t = s, 
like that of the parton distribution function fi(x,t) with t = s (see Section 35 of U), is 
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governed by the DGLAP equation M 

t^D^t) = Y / J^^PAz,a s )D,(x/z,t) . (5) 

In analogy to DIS, in some cases an evolution equation for the fragmentation function 
F itself (Eq. (|j)) can be derived from Eq. [Q. Notice that the splitting function is 
now Pji rather than Pij since here Dj represents the fragmentation of the final parton, 
whereas fj in Eqs. (35.7) - (37.10) of [|TJ represented the distribution of the initial parton. 
The splitting functions again have perturbative expansions of the form 



P Jt (z,a s ) = P^\z) + ^l\z) + --- (6) 



where the lowest-order functions P^ (z) are the same as those in deep inelastic scattering 



but the higher-order terms |f§[] are different. The effect of evolution is, however, the 
same in both cases: as the scale increases, one observes a scaling violation in which the x 
distribution is shifted towards lower values. This can be seen from Fig. [l]. 

As in the case of the parton distribution functions, the most common strategy for solving 
the evolution equations (|5p is to take moments (Mellin transforms) with respect to x: 

D(j,t)= f 1 dx x*- 1 D(x,t) , (7) 
Jo 

the inverse Mellin transformation being 

D(x,t) = J c dj x~ j D(j,t) , (8) 

where the integration contour C in the complex j plane is parallel to the imaginary axis 
and to the right of all singularities of the integrand. Again we can consider fragmentation 
function combinations which are non-singlet (in flavour space) such as Dy = D q . — D q . 
or D qi — D qj . In these combinations the mixing with the flavour singlet gluons drops out 
and for a fixed value of a s the solution is simply 

D v (j,t)=D v (j,t )(-) , (9) 



where the quantity ^ qq is the Mellin transform of a s P qq /2n. The resulting dependence 
on t violates scaling rules based on naive dimensional analysis, and therefore the power 
7 ?g is called an anomalous dimension. For a running coupling a s (t), however, the scaling 
violation is no longer power-behaved. Inserting the lowest-order term in Eq. (9.5a) of M 
for the running coupling, we find that the solution varies like a power of hit: 

D v (j,t) = 

d qq (j) = -,-,,-r ---^-^7 • (10) 




lr rhere are misprints in the formulae in the published article. The correct expressions can be found 
in the preprint version or in M. 
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Figure 1: The e + e~ fragmentation function for all charged particles is shown (a) for 
different cm. energies, \/s, versus x and (b) for various ranges of x versus \fs. For the 
purpose of plotting (a), a scale factor c(y/s) = 10* was applied to the distributions where 
% is ranging from i = (y/s = 12 GeV) to % = 12 (\/s = 189 GeV). Data are compiled 
from references @, g g g 0, 0, |T|, [14], p], H, 0, |T|, ||. 
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where Cp = 4/3, Ca = 3, and rif is the number of quark flavours as in Section 9 of [T|]. 
For the singlet fragmentation function 

A? = £(A, + £y (ii) 

i 

we have mixing with the fragmentation of the gluon and the evolution equation becomes 
a matrix relation as in the deep inelastic scattering case: 



\ = ( 7« 2n n 9 1 
9 J \ ^99 ^gg / \ ^g 

The anomalous dimension matrix in this equation has two real eigenvalues, 





(12) 



7± 



9 hgg + Iqq ± \j ilgg ~ 7 W ) 2 + 8ri /7 OT 7«,] • (13) 



Expressing D$ and D g as linear combinations of the corresponding eigenvectors D + and 
£)_, we find that they evolve as superpositions of terms of the form ( |i~0| ) with 7+ and 7_ 
in the place of j qq . 

At small x, corresponding to j — > 1, the g ^ g anomalous dimension becomes dominant 
and we find 7+ — > 7 g9 — > 00, 7_ — > 7 99 — > 0. This region requires special treatment, which 
will be presented in Section ^. 

There are several complications in the experimental study of scaling violation in jet frag- 
mentation functions f|. First, the energy dependence of the electroweak couplings gi(s) 
that enter into Eq. @ is especially strong in the energy region presently under study 
(y/s = 20 — 200 GeV). In particular, the 6-quark contribution varies by more than a 
factor of 2 in this range. The fragmentation of the b quark into charged hadrons, includ- 
ing the weak decay products of the 6-flavoured hadron, is expected to be substantially 
softer than that of the other quarks, so its varying contribution can give rise to a 'fake' 
scaling violation. A smaller, partially compensating effect is expected in charm fragmen- 
tation. These effects can be eliminated by extracting the b and c fragmentation functions 
from tagged heavy quark events. Fig. |2| shows the flavour-dependent e + e~ fragmentation 
functions for quarks determined at y/s = 91 GeV. 

Secondly, one requires the gluon fragmentation function D g (x,s) in addition to those of 
the quarks. Although the gluon does not couple directly to the electroweak current, it 
contributes in higher order, and mixes with the quarks through evolution. Its fragmenta- 
tion can be studied in tagged heavy-quark three-jet (QQg) events, or via the longitudinal 
fragmentation function (to be discussed below). Results for both methods are shown in 
Fig. |. 

A final complication is that power corrections to fragmentation functions, of the form 
f(x)/Q p , are not so well understood as those in deep inelastic scattering. Theoretical 
arguments [^, ^jj suggest that hadronization can lead to 1/s corrections. Therefore, 
possible contributions of this form should be included in the parametrization when fitting 
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Figure 2: Comparison of the charged-particle and the flavour-dependent e + e~ fragmen- 
tation functions obtained at y/s = 91 GeV. The data are shown (a) for the inclusive, light 
(up, down, strange) quarks, charm quark, bottom quark, and the gluon versus x, and (b) 
in various ranges of x. For the purpose of plotting (a), a scale factor c(flavour) = 10* was 
applied to the distributions where % is ranging from i = (Gluon) to i = 4 (all flavours). 



Data are compiled from references [12, [L3l 14, 15, 16 



the scaling violation. More conservatively, one could also try to fit 1/Q corrections to 
check for their presence. 

Quantitative results of studies of scaling violation in e + e~ fragmentation are reported in 
311. The values of a s obtained are consistent with the world average 



refs. H, H, [ 
in Section 9 of 111 



3 Average Multiplicities 

The average number of hadrons of type h in the fragmentation of a parton of type i at 
scale t, (Tihit))^ is just the integral of the fragmentation function, which is the j ' — 1 
moment in the notation of Eq. (^): 

(n h (t)) t = [ l dxD?(x,t) = D?(l,t). (14) 

If we try to compute the t dependence of this quantity, we immediately encounter the 
problem that the lowest-order expressions for the anomalous dimensions 7 OT and 7 99 in 
Eq. ( |T2| ) are divergent at j — 1. The reason is that for j < 1 the moments of the splitting 
functions in Eq. @ are dominated by the region of small z, where P g i{z) has a divergence 
associated with soft gluon emission. 
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various experiments is shown [O, O, O, 



The mean multiplicity of charged particles measured in e + e 

[I7L pi pi M m H, M, [37L 138 



annihilation by 
91 BOiP, 



|2, [f^, KB[ . The measurements include contributions from Kg and A decays. 

Overlaid is the prediction Eq. (|16|) using a s (M|) = 0.118. 



In fact, however, we can still solve the evolution equation for the average multiplicity 
provided we take into account the suppression of soft gluon emission due to coherence 
|32|, [33]. The leading effect of coherence is that the scale on the right-hand side of the 
DGLAP equation Eq. (||) is reduced by a factor of z 2 : 



" l dz a 



(15) 



This change is not important for most values of x but it is crucial at small x. The 
anomalous dimensions are now found to be finite at j = 1 and Eq. (0) gives to next-to- 



leading order p4| : 



(n(s)) = a ■ exp 




+ (i + H) lnQs(s) 



+ c 



(16) 



where a and c are constants, and f3 Q is defined in Eq. (9.4b) of [|I|]. The resulting prediction, 
shown by the curve in Fig. [| after fitting the parameters a and c, is in very good agreement 
with experiment. In this comparison the scale parameter A mu i t in a s (s) was held fixed. 
A mu it is not necessarily equal to A^g, because the renormalization scheme dependence 
of (n(s)) only appears at next-to-next-to- leading order, and corrections to Eq. flTtjD of 
this order have not yet been calculated. In fact, however, the value used in Fig. 0, 
A mu it = 226 MeV, corresponding to a s (M|) = 0.118, is close to A^g, indicating that 
further higher-order corrections should be small. Higher order corrections to Eq. (|T6|) 

1,0- 



have been considered in 48 
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4 Small- £ Fragmentation 



The behaviour of D(j, s), Eq. ([?]), near j — 1 determines the form of small- x fragmentation 
functions via the inverse Mellin transformation (§). Keeping the first three terms in a 
Taylor expansion of the anomalous dimension 7 99 around j = 1 gives a simple Gaussian 
function of j which transforms into a Gaussian in the variable £ = ln(l/x): 

1 



xD(x, s) oc exp 



2a 2 



where the peak position is 



1 



Aba s (s) 

and the width of the distribution of £ is 



In 



a 




2tt 



oc 



s 

A2 



(17) 



(18) 



(19) 



C A a|(s), 

Again, one can compute next-to-leading corrections to these predictions. In the method 



of pniifl , the corrections are included in an analytical form known as the 'modified leading 
logarithmic approximation' (MLLA). Alternatively |5l| they can be used to compute the 
higher moment corrections to the Gaussian form (|T7|). Fig. ^ shows the £ distribution at 
various cm. energies. 



The predicted energy dependence (|18|) of the peak in the £ distribution is a striking 
illustration of soft gluon coherence, which is the origin of the suppression of hadron 
production at small x. Of course, a decrease at very small x is expected on purely 
kinematical grounds, but this would occur at particle energies proportional to their masses, 
i.e. at x oc m/y/s and hence £ ~ | In s. Thus if the suppression were purely kinematic the 
peak position £ p would vary twice as rapidly with energy, which is ruled out by the data 
(see Fig. |5|). 



5 Large-x Fragmentation 

At large values of the energy fraction x, there are enhancements in the coefficient functions 
Ci in (|) and the splitting functions Pji in Eq. (|5|). These are associated with the emission 
of soft and/or collinear gluons, which can lead to factors of a" ln m (l — x) /(1 — x) with m < 
In — 1 in the nth order of perturbation theory. It turns out that in the conventional MS 
factorization scheme the largest terms, with m > n, all occur in the coefficient functions. 
After the Mellin transformation (|7|), they are resummed by the following expression pD , |61[1 

lnC g (j,s) = - ['dz^—^ I f ^A[a s (t)} + B[a s ((l - z)s}\ (20) 

JO 1 — Z [J(l-z)s t J 

where A(a s ) and B(a s ) have perturbative expansions 

M«s) = E f ^y^ (n) , Bias) = £ (^) n B^ (21) 



n=l ^ 71 / n=l V 71 
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Figure 5: Evolution of the peak position, £ p , of the £ distribution with the cm. energy y/s. 
The MLLA QCD prediction (solid) and the expectation without gluon coherence (dashed) 
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were fitted to the data 
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with 



c F --, 



c F 

3 



67 
36 



71 

12 



67 7T 2 10 



9 



27 



(22) 



C f 



6 Longitudinal Fragmentation 



In the process e + e~ — > — ► ftjf , the joint distribution in the energy fraction x and the 
angle between the observed hadron h and the incoming electron beam has the general 
form 

1 d 2 a 3 3 . o „ „ , s 3 



0V 



~(1 + coar 6) F T (x) + -sia 2 9F L (x) + -cos9F A (x) 



(23) 



_ _dxd cos 9 8 y ' ' 4 " K ' 4 

where Ft, Ft and Fa are respectively the transverse, longitudinal and asymmetric frag- 
mentation functions. All these functions also depend on the cm. energy y/s. Eq. (|23|) 
is the most general form of the inclusive single particle production from the decay of a 
massive vector boson [|J. As their names imply, Ft and Ft represent the contributions 
from virtual bosons polarized transversely or longitudinally with respect to the direction 
of motion of the hadron h. Fa is a parity-violating contribution which comes from the 
interference between the vector and axial vector contributions. Integrating over all angles, 
we obtain the total fragmentation function, F = Ft + Ft- Each of these functions can be 
represented as a convolution of the parton fragmentation functions Di with appropriate 
coefficient functions Cj ,L ' A as in Eq. (^). This representation works in the high energy 
limit. Furthermore, when x ■ \/s is of the order of a few hundred MeV, the p± smear- 
ing of the parton can no longer be neglected, and the fragmentation function formalism 
no longer account correctly for the separation of Ft, Ft, and Fa- The transverse and 
longitudinal coefficient functions are 



C h q {z) 



5{l-z) + 0{a s ) 
0(a s ) 



0(a 



which implies that 

F L (x) 



1 dz 



2tt Jx z 



F t (z) 



1 l Dg( Z ) 



0(ai 



(24) 



(25) 



Thus the gluon fragmentation function D g can be extracted from measurements of Ft 
and Ft- The next-to-leading order corrections to Eq. ([25|) are also known |)4], p5]| . In 
Fig. |6| F T , Ft, and Fa measured at yfs = 91 GeV are shown. The gluon fragmentation 
function derived from Ft and Ft in [24 is shown in Fig. 0(a). 
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Figure 6: Transverse {Ft), longitudinal (Fl), and asymmetric {Fa) fragmentation func- 



tions are shown [13[ |24], |29|. Data points with relative errors greater than 100% are not 
shown. 
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Summed over all particle types, the total fragmentation function satisfies the energy sum 
rule (^|), which we may write as 

1 r l 

alxxF(x) = 1 . (26) 



Similarly the integrals 

1 f 1 dxxF T . L (x) = ^ (27) 

2 JO cr tot 

give the transverse and longitudinal fractions of the total cross section. The perturbative 
prediction is |64 

a L a s /601 6 37 \ /«s\ 2 . n( ^ 



^ + (13.583 - 1.028^)^^ (28 



where £(3) = 1.202. . .. Comparing with Eq. (9.17) of Q, we see that the whole of the 
0(a s ) correction to a tot comes from the longitudinal part, while the 0(af) correction 
receives both longitudinal and transverse contributions. 



7 Gluon fragmentation 

As we saw in the previous Section, the gluon fragmentation function can be extracted 
from the longitudinal fragmentation function using Eq. ( |25"D . It can also be deduced from 
the fragmentation of three-jet events in which the gluon jet is identified, for example by 
tagging the other two jets with heavy quark decays. The trouble with the latter method is 
that the relevant coefficient function has not been computed to next-to-leading order, and 
so the scale and scheme dependence of the extracted fragmentation function is ambiguous. 
The experimentally measured gluon fragmentation functions are shown in Fig. 0(a). 



8 Fragmentation models 



Although the scaling violation can be calculated perturbatively, the actual form of the 
parton fragmentation functions is non-perturbative. Perturbative evolution gives rise to 
a shower of quarks and gluons (partons). Phenomenological schemes are then used to 
model the carry-over of parton momenta and flavour to the hadrons. Two of the very 
popular models are the string fragmentation [^, |67|, implemented in the JETSET ||| 



and UCLA [^] Monte Carlo event generation programs, and the cluster fragmentation of 



the HERWIG Monte Carlo event generator |70 |. 
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8.1 String fragmentation 



The string- fragmentation scheme considers the colour field between the partons, i.e., 
quarks and gluons, to be the fragmenting entity rather than the partons themselves. 
The string can be viewed as a colour flux tube formed by gluon self-interaction as two 
coloured partons move apart. Energetic gluon emission is regarded as energy-momentum 
carrying "kinks" on the string. When the energy stored in the string is sufficient, a qq 
pair may be created from the vacuum. Thus the string breaks up repeatedly into colour 
singlet systems as long as the invariant mass of the string pieces exceeds the on-shell mass 
of a hadron. The qq pairs are created according to the probability of a tunneling process 
exp(— 7rm^ ± / k) which depends on the transverse mass squared rn^ ± = m 2 q +Pg ± and the 
string tension k « 1 GeV/fm. The transverse momentum p Qy ± is locally compensated 
between quark and antiquark. Due to the dependence on the parton mass m q and/or the 
hadron mass m/,, the production of strange and, in particular, heavy-quark hadrons is 
suppressed. The light-cone momentum fraction z = (E +p\\)f l /(E +p) q , where p\\ is the 
momentum of the formed hadron h along the direction of the quark q, is given by the 
string-fragmentation function 

/(*0~i(l-z) a exp(-^) (29) 

where a and b are free parameters. These parameters need to be adjusted to bring the 
fragmentation into accordance with measured data, e.g. a = 0.11 and b = 0.52 GeV -2 as 
determined in JFI] (for an overview see |ff2|). 



8.2 Cluster fragmentation 

Assuming a local compensation of colour based on the pre- confinement property of pertur- 
bative QCD [fT3"| , the remaining gluons at the end of the parton shower evolution are split 
non-perturbatively into quark- antiquark pairs. Colour singlet clusters of typical mass of a 
couple of GeV are then formed from quark and antiquark of colour-connected splittings. 
These clusters decay directly into two hadrons unless they are either too heavy (relative 
to an adjustable parameter CLMAX, default value 3.35 GeV), when they decay into two 
clusters, or too light, in which case a cluster decays into a single hadron, requiring a small 
rearrangement of energy and momentum with neighbouring clusters. The decay of a clus- 
ter into two hadrons is assumed to be isotropic in the rest frame of the cluster except if 
a perturbative-formed quark is involved. A decay channel is chosen based on the phase- 
space probability, the density of states, and the spin degeneracy of the hadrons. Cluster 
fragmentation has a compact description with few parameters, due to the phase-space 
dominance in the hadron formation. 
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9 Experimental studies 



A great wealth of measurements of e + e~ fragmentation into identified particles exists. 
A collection of references, where data on the fragmentation into identified particles can 
be found, is given for Tab. 37.1 of |]lj. As representatives of all the data, Fig. [7| shows 
fragmentation functions as the scaled momentum spectra of charged particles at several 
cm. energies, y/s. In Fig. |8| x p spectra at y/s = 91 GeV are shown for all charged particles, 
for several identified charged and neutral particles. Heavy flavour particles are dealt with 



separately in Sect. |10 



The measured fragmentation functions are solutions to the DGLAP equation ([5]) but need 
to be parametrized at some initial scale to (usually 2 GeV 2 for light quarks and gluons). 
A general parametrization is |7£ 



D^ h (x,t Q ) = Nx a (l-x) f> ( 



1 + 1 

x 



(30) 



81 



83 



where the normalization N, and the parameters a, /3, and 7 in general depend on the 
energy scale to and also on the type of the parton, p, and the hadron, h. Frequently 
the term involving 7 is left out |65], |8 



In the above quoted references the 
(]) are tabulated for many different combinations of partons and 
The parameters were obtained by fitting data for many different 

5 - 200 GeV). 



parameters of Eq. 
hadrons in p — > h. 

hadron types over a vast range of cm. energies (i/i 



10 Heavy quark fragmentation 



It was recognized very early []84j that a heavy flavoured meson should retain a large 



fraction of the momentum of the primordial heavy quark, and therefore its fragmentation 
function should be much harder than that of a light hadron. In the limit of a very heavy 
quark, one expects the fragmentation function for a heavy quark to go into any heavy 
hadron to be peaked near 1. 

When the heavy quark is produced at a momentum much larger than its mass, one expects 
important perturbative effects, enhanced by powers of the logarithm of the transverse mo- 
mentum over the heavy quark mass, to intervene and modify the shape of the fragmenta- 
tion function. In leading logarithmic order (i.e., including all powers of a s \ogmq/pT) the 
total (i.e., summed over all hadron types) perturbative fragmentation function is simply 
obtained by solving the leading evolution equation for fragmentation functions, Eq. (|5J), 
with the initial condition at a scale fi 2 = mq given by 

D Q (z,m 2 Q ) = 5(l-z), D q = 0, Dg = , D Q = , D g = (31) 

where the notation D^z) stands now for the probability to produce a heavy quark Q from 
parton i with a fraction z of the parton momentum. If the scale fi 2 is not too large, and 
one looks at relatively large values of z, one can assume that the non- singlet evolution 
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embodies most of the physics. Following this assumption, for example, the average value 
of z at a scale /i 2 is easily obtained by solving the evolution equation in the moment 
representation (cf. Eq. (|T0|)) 

SC F 
UC A -2n f 

(32) 

Several extensions of the leading logarithmic result have appeared in the literature: 



\a s {m 2 Q ) 



Next-to-leading-log (NLL) order results for the perturbative heavy quark fragmen- 
tation function have been obtained in 



At large z, phase space for gluon radiation is suppressed. This exposes large per- 
turbative corrections due to the incomplete cancellation of real gluon radiation and 
virtual gluon exchange (Sudakov effects), which should be resummed in order to 
get accurate results. A leading-log (LL) resummation formula has been obtained in 
]B5[ , |86[| . Next-to-leading-log resummation has been performed in 



Fixed order calculations of the fragmentation function at order a 2 in e + e annihila- 
tion have appeared in [37]. This result does not include terms of order (a s log s/m 2 ^ k 



and a s (a s log s/m 2 ) k , but it does include correctly all terms up to the order a 2 , in- 
cluding terms without any logarithmic enhancements. 



Inclusion of non-perturbative effects in the calculation of the heavy quark fragmentation 
function is done in practice by convolving the perturbative result with a phenomenological 
non-perturbative form. Among the most popular parametrizations we have the following: 



Peterson et al. [|S 
Kartvelishvili et al. fBO 



D np (z) oc 
D np (z) oc 



Collins&Spiller [0: D np (z) oc 



z \ z 1 — z 
z a (l-z), 



1 



z 



z)e c 



[l + z 2 



1 

z 



1 - z 



(33) 
(34) 



(35) 



where e, a, and ec are non-perturbative parameters, depending upon the heavy hadron 
considered. In general, the non-perturbative parameters do not have an absolute meaning. 
They are fitted together with some model of hard radiation, which can be either a shower 
Monte Carlo, a leading-log or NLL calculation (which may or may not include Sudakov 
resummation), or a fixed order calculation. In P7|, for example, the e parameter for 



charm and bottom production is fitted from the measured distributions of refs. Wr \)2 



for charm, and of for bottom. If the leading-logarithmic approximation (LLA) is used 
for the perturbative part, one finds e c ~ 0.05 and eb ~ 0.006; if a second order calculation 
is used one finds e c ~ 0.035 and e;, ~ 0.0033; if a NLLO calculation is used instead one 
finds e c ~ 0.022 and e& ~ 0.0023. The larger values found in the LL approximation are 
consistent with what is obtained in the context of parton shower models |94j| , as expected. 
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Figure 9: Efficiency-corrected inclusive cross-section measurements for the production of 
D° and D* + in e + e~ measurements at y/s ~ 10 GeV. The variable x p approximates the 
Peterson variable z, but is not identical to it. 



The e parameter for charm and bottom scales roughly with the inverse square of the heavy 



flavour mass. This behaviour can be justified by several arguments [S4 , p5| , |96| . It can 
be used to relate the non-perturbative parts of the fragmentation functions of charm and 
bottom quarks |87], |97|, |98| . 



The bulk of the available fragmentation function data on charmed mesons (excluding 
J/^(1S)) is from measurements at \fs = 10 GeV.Q Shown in Fig. |9| are the efficiency- 
corrected (but not branching ratio corrected) CLEO |)9| and ARGUS |)2| inclusive cross- 
sections (s • Bda/dxp in units of GeV 2 nb, with x p = p/p max ) for the production of pseu- 
doscalar D° and vector D* + in e + e~ annihilation at « 10 GeV. For the D°, B rep- 
resents the product branching fraction: D* + — > D 7T + , D° — > K~n + . These inclusive 
spectra have not been corrected for cascades from higher states, nor for radiative effects. 
Note that since the momentum spectra are sensitive to QED and QCD radiative correc- 
tions, charm spectra at y/s = 10 GeV cannot be compared directly with spectra at higher 
cm. energies, and must be appropriately evolved. 

Tuning e of the function (|33D in the JETSET Monte Carlo generator pHf using the pa- 



rameter set of [71] and including radiative corrections to describe the combined CLEO 
and ARGUS D° and D* + data gives e c = 0.043 ± 0.004 (x 2 /d.o.f.= 50.5/41); this is in- 
dicated in the solid curves. Measurements of the fragmentation functions for a variety of 
particles has allowed comparisons between mesons and baryons, and particles of differ- 
ent spin structure, as shown in Table [L| The e values listed in this table were obtained 
from directly fitting the Peterson function (|33|) to the measured differential cross-section 
cr • Bda/dXp. 

2 This part on charm quark fragmentation including Fig. ^ and Table [j], which have been updated, are 
taken from the review of D. Besson contributed to C. Caso et al., Eur. Phys. J. C3, 1 (1998). 
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Table 1: The Peterson momentum hardness parameter e as obtained from fits of Eq. (|33|) 
to s ■ Bda/dxp of e + e~ — ► particle + X measurements at ^/s ~ 10 GeV. 

Particle 



D 

D+ 
D*+ 
D s 
D* 



0.260 ±0.024 

0.156 ±0.022 

0.198 ±0.022 

0.10 ±0.02 

0.057 ±0.008 



\Y.)\ \\{\{\ 



sra 
Bra 



D 1 (2420) u 

L>*(2460)° 

Di(2420)+ 

L>*(2460) + 

D sl (2536) 4 

^j(2573)- 



1 0.015 ±0.004 

1 0.039 ±0.013 

1 0.013 ±0.006 

1 0.023 ±0.009 

1 0.018 ±0.008 

1 0.027^ 



K).043 
-0.015 
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loc| 
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10C| 


102. 


loq 


102, 


locj 


103, 


10C| 


103. 


10G| 


104, 


100| 



in 



' — c 





E c (2455) 
S c (2520) 
S C (2645) H 
H c (2645) c 













0.267 ±0.038 
0.23 ±0.05 

0.28 ±0.05 

so +0 - 10 

0.24±g" 22 



106, 107 



[ 1081 TO 



E3 



[mj, gg 



0.22 



A c (2593)- 
A c (2625)" 
S c (2815) 



10 
0.15 
-0.08 



1 0.058 ±0.022 
1 0.053 ±0.014 
1 0.07^ 



EH 
[rjjl 
E3 



116j |H7 



K).03 
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EH 
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Figure 10: Measured e + e fragmentation function of b quarks into B hadrons ||121|, |122 



123 . 



We note from Table |l]that the mass dependence of e is less marked than the dependence on 
the orbital angular momentum structure of the charmed hadron being measured. Orbitally 
excited L = 1 charmed hadrons (Dj, D s j, and A c (2593), A c (2625)) show consistently 
harder spectra (i.e., smaller values of e) than the L = ground states, whereas the data 
for the ground state charmed baryons A c and S c show agreement with the lighter (by 
~ 400-600 MeV) ground-state D and D s charmed mesons. To some extent, the harder 
spectra of L = 1 hadrons can be attributed to the fact that all the L = 1 charmed hadrons 
will eventually decay into L = hadrons. 

Experimental studies of the fragmentation function for b quarks have been performed at 
LEP and SLD 



3 120, 121, 122 



The results are shown in Fig. [T(J. B hadrons are com- 
monly identified by their semileptonic decays or by their decay vertex reconstructed from 
the charged particles emerging from the B decay. The most recent results of refs. 



123 



and ||122j1 are in agreement. Both studies fit the B spectrum using a Monte Carlo shower 
model supplemented with non-perturbative fragmentation functions to fit their data. 

The experiments measure primarily the spectrum of B meson. This defines a fragmenta- 
tion function which includes the effect of the decay of higher mass excitations, like the B* 
and B**. There is an ambiguity in the definition of the fragmentation function, which has 
to do with what is considered to be a final state hadron. One may prefer to distinguish 
B meson produced directly from those arising from decays of B* and B**, and define 



(36) 



where the notation stands for directly produced. The fragmentation function which 
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Table 2: Fraction of events containing g — > cc and g — > bb subprocesses in Z decays, 
as measured by the various collaborations, compared with theoretical predictions. The 
central/lower /upper values for the theoretical predictions are obtained with m c = (1.5 ± 
0.3) and m b = (4.75 ± 0.25) GeV. 



Ug^cc (%) W^fcg (%)' 



OPAL 
ALEPH 
DELPHI 
SLD 


128 
129 


3.20 ±0.21 ±0.38 
2.65 ±0.74 ±0.51 


PIT 

m 

132 


0.277 ±0.042 ±0.057 
0.21 ±0.11 ±0.09 
0.307 ±0.071 ±0.066 


Theory 126 





A^= 150 MeV 1.35™ 0.20 ± 0.02 



Ai^ = 300MeV 1.85™ 0.26 ± 0.03 



is directly measured from the B mesons, irrespective of their origin, is instead given by 

D B (x) = Df{x) ± J F B ^ B (y)D^(z)5(x - yz)dydz 

+ J F B ^ B (y)D B i l(z)8(x - yz)dydz . (37) 

where (for example) F B **^ B (y) stands for the probability for a fast B** to decay to B with 
a fraction y of its energy. Thus, in order to extract -D^b*|b« from D B a, correction factor 
should be computed, which gives a harder fragmentation function. Thus, for example, 
ref. pi obtains the value (x B ) = 0.714 ± 0.005(stat.) ± 0.007(syst.) ± 0.002(model); 



postulating a B* fraction of 0.75, one gets {x B \ B *) = 0.718; postulating a B** fraction of 
0.25 yields (x BlB **) = 0.728. 

Besides degrading the fragmentation function by gluon radiation, QCD evolution can also 
generate soft heavy quarks, increasing in the small x region as s increases. This effect has 
been studied both theoretically and experimentally. One important issue is to understand 
how often bb or cc pairs are produced indirectly, via a gluon splitting mechanism. Several 
theoretical studies are available on this topic | |124| , |125| , |126| , |127| ]. Experimental studies 



on charm production via gluon splitting have been presented in refs. [|128] , |129| , and 



measurements of g — > bb have been given in ||130| , |131| , |132|| . The reported values are given 



in Table |2|. In ref. ||126|| an explicit calculation of these quantities has been performed. 
Using these results charm and bottom multiplicities for different values of the masses and 
of A^- were computed in ||133|| . They are reported in Table |2|. As can be seen, charm 
measurements are somewhat in excess of the predictions. The averaged experimental 
result for charm, (3.10 ± 0.34)%, is 1 - 2 standard deviations outside the range of the 
theoretical prediction, preferring lower values of the quark mass and/or a larger value of 
A^-. However, higher-order corrections may well be substantial at the charm quark mass 
scale. Better agreement is achieved for bottom. 

As reported in ref. [|126|1 , Monte Carlo models are in qualitative agreement with these 
results, although the spread of the values they obtain is somewhat larger than the theo- 
retical error estimated by the direct calculation. In particular, for charm one finds that 
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while HERWIG 
ARIADNE " 



(J and JETSET agree quite well with the theoretical calculation, 



134 1 is higher by roughly a factor of 2, and thus is in better agreement with 



data. For bottom, agreement between theory, models and data is adequate. For a detailed 
discussion see ref. 



135] 



The discrepancy with the charm prediction may be due to experimental cuts forcing the 
final state configuration to be more 3-jet like, which increases the charm multiplicity. 



Calculations that take this possibility into account are given in [127] 
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